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ABSTRACT 

In this research a Collaborative Optimization (CO) approach 
for multidisciplinary systems design is used to develop a 
decision based design framework for non-deterministic 
optimization. To date CO strategies have been developed for use 
in application to deterministic systems design problems. In this 
research the decision based design (DBD) framework proposed 
by Hazelrigg (1996a, 1998) is modified for use in a 

collaborative optimization framework. The Hazelrigg 
framework as originally proposed provides a single level 
optimization strategy that combines engineering decisions with 
business decisions in a single level optimization. By 
transforming this framework for use in collaborative 
optimization one can decompose the business and engineering 
decision making processes. In the new multilevel framework of 
Decision Based Collaborative Optimization (DBCO) the 
business decisions arc made at the system level. These business 
decisions result in a set of engineering performance targets that 
disciplinary engineering design teams seek to satisfy as part of 
subspace optimizations. The Decision Based Collaborative 
Optimization framework more accurately models the existing 
relationship between business and engineering in 

multidisciplinary systems design. 


NOMENCLATURE 

C F fixed cost 

C r total cost 

C unit cost 

II 

C v variable cost 

£'(?) expected value 

F objective function in optimization problem 

NR net revenue 

p 

price 

PNR present value of net revenue 

a vector of attributes 

,;v demand coefficient for attribute / 

demand coefficient for price 

price r 

d vector of compatibility (or discrepancy) constraint 

d i compatibility constraint in discipline / 

g vector of constraints 

g. vector of constraints in discipline i 

u number of cost-related variables 

n number of attributes 

a 



number of subspaccs (disciplines) 
vector of parameters 

variable cost parameter associated with cost-related 

variable i 
demand 

baseline demand 
utility 

vector of design variables 
vector of auxiliary design variables 
vector of auxiliary design variables corresponding to 

vector of local design variables for discipline i 
vector of shared design variables 
vector of shared design variables used in discipline i 
vector of subspace design variables in discipline / 

vector of optimal target values in discipline i 
vector of system level design variables 

vector of targets (i.e., system level design variables) 
sent down to discipline / 
vector of state variables 

vector of state variables calculated as output from 
discipline / 

vector of coupling variables calculated as output from 
discipline i and used as input to discipline j 
system level target variable 
variable at the optimum design 

INTRODUCTION 

Increasing attention has been paid to the notion that 
engineering design is a decision-making process (Chen, ct. 
al, 1 998, 2000, Kim, cl. al, 2000, Azarin, 2000). This notion is 
consistent with the definition of decision as a choice from 
among a set of options and as an irrevocable allocation of 
resources. The approach of decision-based design (DBD) is 
built upon this notion. Rooted from more than two hundred 
years of research in the field of decision science, economics, 
operations research and other disciplines, decision-based design 
(DBD) provides a rigorous foundation for design, which enables 
engineers to identify the best trade-off and focus on where the 
payoffs are greatest. 

Engineering design involves the generation of design 
alternatives or options and the selection of the best one. Since 
the number of possible design options is practically infinite for 
most products, human judgement is needed to decide which 
options to include in the consideration of alternative designs and 
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which to neglect. Moreover an appropriate value measure has to 
be determined in order to compare and rank order design 
options. It is impossible to know exactly how a particular design 
alternative will perform before it is built. However the product 
cannot be built until it is selected. Evidently engineers have to 
make their selection a priori , without full knowledge of the 
consequence of this certain selection. Thus design is always a 
matter of normative decision making under uncertainty and risk. 

Optimization techniques have been widely applied to select 
‘he preferred design options from the set of alternatives taken 
into consideration without having to explicitly evaluate all 
possible alternatives in the set. There exists a close relationship 
between decision making and optimization. In general, decision 
making has three elements: generation or identification of 
options, assignment of expectations on each option and the 
application of preferences to determine the preferred choice. An 
optimization problem involves the maximization or 
minimization of an objective function or functional F(x) in the 
feasible region of design variables jc. A careful comparison 
between decision making and optimization will reveal that 
options, expectations and preferences are all present in 
optimization. The option space is equivalent to the set ot 
permissible values of x in the feasible region. The expectation 
of any given x is assigned by F(x)> and the preference is stated 
That more is better (maximization) or less is better 
(minimization). Thus optimization can be used to capture the 
properties of decision-making. This recognition allows the 
application of rigorous decision theory to the case of 
optimization. 

DECISION-BASED DESIGN (DBD) FRAMEWORK 

Application of decision-based design within an optimization 
domain, requires practitioners to formulate valid objective 
functions for proper decision-making. Most of the research in 
the field of optimization has been focused on the solution ot the 
optimization problem such as development, improvement and 
implementation of search methods to locate the optimum, while 
little attention has been paid to optimization problem 
formulation. In fact, the issue of problem formulation is of the 
same significance, if not more, as the issue of finding the 
optimal solution. An solution obtained using any search method 
is no better than the objective function chosen lor the 
optimization. If an irrelevant objective is used, the solution is 
equally irrelevant (Hazclrigg 1996a, 1997). Therefore a primary 
concern in DBD is the development a mathematically sound 
objective function. Recognizing that design is a decision- 
making process, it is imperative to construct an objective 
function which is able to capture the preferences of rational 
decision makers in system design problems involving risk. The 
decision-based design (DBD) framework of Hazclrigg 1996a, 
1998 (Fig. 1) provides a basis for exploring this issue. 
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The Rule of Rational Decisions 

The aforementioned decision-based design framework 
implements the concept of rational decisions. Rational decisions 
follow the rule that the preferred decision is the option whose 
expectation has the highest value. In a normative approach, 
decisions involve options, expectations and values. An 
expectation combines the possible outcomes of an alternative 
with probabilities of occurrence of each possible outcome. Note 
that in general expectations are not equivalent to outcomes. An 
outcome refers to what actually happens after a decision is made 
to select a certain option, while an expectation refers to what is 
expected to happen, based on available knowledge, as the result 
of a decision before the decision is made. In other words 
expectations are associated with what will happen, therefore, 
expectations relate to the future. In the process of engineering 
design, it is practically impossible to predict the future with 
precision and certainty. The outcomes of most options in the 
design option space cannot be determined with certainty prior to 
the decision to select one option. Hence expectations are always 
probabilistic. Engineers are forced to make decisions under 
uncertainty and risk. 

Unlike problem-solving, decision-making cannot be 
conducted in the absence of human values. The purpose of 
values in decision making is to rank order alternatives according 
to the preference of decision makers. In the case of 
optimization, engineers seek a design which maximizes value, 
and the objective function serves as a numerical value function 
to automate the process of rank ordering. By means of an 
objective function, a real scalar is assigned to each design 
alternative in accordance with the decision maker’s preference. 
In this sense an objective function serves as a utility function in 
the context of economics. Note that utilities are determined by 


preferences and so is the objective function. A necessary 
condition for the existence of an objective function in decision- 
based design is that decision makers or design engineers are 
rational individuals whose preferences and indifferences 
between all pairs of outcomes in the design space exists and 
comprises a transitive set. 

Decision Making/Optimization Under Uncertainty 

Due to the nature of engineering design, expectations on 
design alternatives can never be determined with certainty. 
Particular care should be taken in the formulation of an 
objective function when risk is present. It is imperative that the 
objective function (or utility function in the context of 
economics) must be valid under conditions of uncertainty and 
risk. The von Neumann-Morgenstem (vN-M) utility is such a 
aluc measure. Built upon the notion of von Neumann- 
Morgenstem lottery and six rigorous axioms (von Neumann, 
953), the normative framework for decision making under risk 
leads to a simple but profound result — the so-called expected 
utility theorem: “The utility of a lottery is the sum of the utilities 
of all possible outcomes of the lottery weighted by their 
probabilities of occurrence.” It follows that the preferred choice 
Tom among a set of risky options is the option with the highest 
expected utility. In the case of optimization, the axioms of vN- 
M utility should be adhered to and the objective function should 
entail the assignment of a vN-M utility to each design 
alternative under consideration. 

Characteristics of DBD Framework 

The nonnative framework of decision-based design 
implements the concept of rational decisions. It facilitates vN-M 
utility as a measure of value against which design alternatives 
ran be compared and optimal designs sought. 



Figure 1. Decision-Based Design Framework (Hazelrigg, 1996b) 
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The DBD framework may view the objective of systems 
design as one of maximizing profit. Profit, which is also referred 
to as net revenue (NR), is revenue generated by the product less 
all costs generated by the product. Revenue is the sum of 
products sold times their prices, in other words, it can be 
calculated as the product of the demand q and the price P at 
which the product is sold. Costs are the sum of things bought 
multiplied by their prices. Total cost Cj consists of: cost of 
manufacture, and all other life cycle costs such as costs of 
research, design, maintenance and repair and much more. The 
calculation of net revenue (NR) can be summarized as Eq. (1): 

NR = P q-C T ( 1 ) 

To account for the time value of money since revenues and 
costs are spread over periods of time, a better index for profit 
would be the present value of net revenue ( PNR ) which can be 
obtained from properly discounting net revenue (NR) to the 
present and integrating it over time. 

Recognizing most products are designed to make money, a 
valid objective function for optimization or decision making 
under uncertainty and risk can be established according to the 
rule of rational decisions: in the process of decision-based 
design, the optimizer should seek to maximize the expected vN- 
M utility of the profit., or net revenue (NR). 

The axioms of vN-M utility dictates that utilities are not 
arbitrarily selected. Rather, utilities must provide proper rank 
ordering of design alternatives in the case of vN-M lottery. In 
non-deterministic optimization where uncertainty and risk 
cannot be neglected, the decision maker’s risk preference 
towards money must be taken into consideration. In general, the 
risk preference of the decision maker towards money leads to a 
utility of money that has diminishing marginal value (Bernoulli, 
1738). Thus profit or net revenue (NR) itself is not a valid utility 
under such circumstances. 

MULTIDISCIPLINARY ENTERPRISE MODEL 

Design is inherently a multidisciplinary process. Traditionally 
multidisciplinary design has focused on disciplines within the 
field of engineering analysis, such as aerodynamics, solid 
mechanics, kinematics, control, and many others. In the 
approach of decision-based design, engineers arc compelled to 
look at the design from a broader viewpoint. Indeed, decision- 
based design encompasses not only the engineering disciplines 
but also business disciplines including economics, marketing, 
operation research and more. In addition to aiming at improving 
the performance of a design, engineers must be aware of the 
substantial impact of non-engineering disciplines on the goal of 
design. Effective communications between engineers and 
experts in the business field are vital to produce successful 
designs. 


The DBD framework of Hazelrigg (1996, 1998) combines 
both engineering and business performance simulations in a 
single level all-at-once optimization approach. In the current 
research the Hazelrigg framework (Fig. 1) has been 
decomposed into the multidisciplinary enterprise model shown 
in Fig. 2. The decomposed system consists of two major 
elements: the engineering disciplines and the business 

discipline. The work in engineering disciplines focuses on 
predicting the performance of the product for different design 
configurations, as well as satisfying performance targets set in 
the business discipline (i.e., management). The role of the 
business discipline centers on providing targets for performance 
improvements in order to yield higher profit. These two 
organizations are coupled through attributes a , total cost C^and 
demand q. 

Attributes a refer to the performances of a product. These 
performances directly influence customer demand. Examples of 
attributes include speed, acceleration, comfort, quality, 
reliability, safety, etc. Prediction of attributes a can be obtained 
from the engineering disciplines. Such information will be used 
as input into the business discipline to estimate the demand q for 
the product, which will then be used to determine the profit or 
net revenue. In order to compute profit, an estimate of the total 
cost Cj is required. C T is again computed in the engineering 
disciplines. Since total cost is influenced by the number of 
product manufactured and/or sold, estimation of demand q is 
called for to feed back into the engineering disciplines. In this 
research it is assumed that the demand for the product, the 
amount of the product manufactured and the amount of the 
product sold are equal. 

The prediction of product attributes a and total cost C T is 
conducted in the engineering disciplines. The engineering 
disciplines module is partitioned into a module o( product 
performance and a module of manufacturing cost and life cycle 
costs. The set of variables which uniquely defines a specific 
design are referred to as engineering design variables x. 
Engineers usually have control over these variables. The 
estimation of demand q y net revenue NR and, ultimately, the 
expected utility of net revenue, arc managed in the business 
discipline. The price of the product not only directly affects the 
amount of profit or net revenue (sec Eq. 1), it is also an 
important factor dnving the demand q , which is also affected by 
the product attributes a. Mathematically speaking, q can be 
modeled as a function of attributes a and the price, P. 

q = q(a , P) (2) 

Note that price is free to be chosen by the decision maker, 
therefore it should be treated as a design variable. 

In the context of multidisciplinary design, the term “system 
analysis (SA)” is often used to describe the process of 
predicting the performance of an engineering artifact using 
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numerical simulation, which is part of the multidisciplinary 
enterprise model investigated in this research. In fact, 
simulation based design tools can be used to assist in the 
analyses within the entire multidisciplinary system. Under these 
circumstances, the analysis of the enterprise model is no 
different than a “system analysis/' only of a larger scale. The 
research results of multidisciplinary design can be readily 
applied to the multidisciplinary enterprise model. Note that the 
equations of physical laws involved in the system analysis arc 
sometimes referred to as “constraints’ 1 in the field of decision 
science. 

The performance predictions obtained from a system analysis 
(SA) are referred to as states y in the context of 
multidisciplinary design. If we think of the multidisciplinary 
enterprise model (Fig. 2) as one big SA, then the attributes a, 
total cost C T , demand q and net revenue NR are components of 
the system states. In this paper states >> are used to represent the 
engineering performance predictions, calculated in the product 
performance module of the multidisciplinary enterprise model. 
Therefore, not all of the states y are attributes a. Only those 
states which influence product demand q will be included in the 
set of attributes a. In decision-based collaborative optimization 
framework (DBCO), the system level optimization treats 
attributes a as system level design variables x* sys which influence 
demand q and serve as targets for the discipline designers to 
satisfy. 


Note that variability exists in the engineering design variables 
x and all simulation-based design tools employed in each 
discipline have prediction errors associated with them. 

SYSTEM OPTIMIZATION 

The multidisciplinary enterprise model discussed above 
oclongs to the type of so-called non-hierarchic systems (or 
ooupled/networked systems). Such systems, usually of fairly 
iarge scale, arc characterized by large numbers of design 
variables x and parameters /?, large numbers of requirements or 
constraints#, and a high level of coupling between participating 
disciplines which are intrinsically linked to one another. 
Coupling results from the information exchange that are 
required within the system analysis. This occurs where the 
output of one discipline is required as input to another discipline 
and vice versa. Some of the mathematical terms, defined by 
Tailing and Sobieszczanski-Sobieski (1994), are briefly 
reviewed below and will be used hereinafter. 

In a coupled system the design variables x can be 
decomposed into the set of shared variables x s h and n S5 sets of 
local variables A*,-, where i ranges from 1 to n SS9 and n ss is the 
iotal number of subspaces (or disciplines). The set of shared 
variables x sfl contains design variables that are needed by more 
ihan one discipline. The set of local variables x, includes design 
variables associated with discipline i only. The set of x^ and the 
>cts of .v, are mutually exclusive subsets of the set of design 


^_.Varia billy Ax 

Engineering Design 
Variables x ^ 



Figure 2. Multidisciplinar y Enterprise Model 
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variables x. The term (x sh )j is used to represent the shared 
design variables that are needed in discipline 

The set of output from discipline i is denoted by y t . The set of 
system states^ is composed of all the j/s in the n ss disciplines. 
The term y t j is introduced to represent the output from discipline 
/' which is also used as input in another discipline j. Note that 
are the subset of coupling variables in the set of 

The vector contains the design constraints associated with 
discipline /. It is assumed that no constraint will stretch over 
more than one discipline. It is also assumed that each inequality 
constraint has been formulated such that zero is the permissible 
value, and the constraint is satisfied when it is greater than or 
equal to zero. In conventional design optimization, constraints 
can be formulated to guard against system failure or to restrict 
the design search to preferred regions of the design space. This 
second class of constraint, related to design preference is not 
used in the decision-based collaborative optimization 
framework developed in this research. In this research preferred 
regions of the design space are imposed implicitly through the 
demand function. 

Collaborative Optimization (CO) 

The Collaborative Optimization (CO) strategy was first 
proposed by Kroo, et al. (1994) and Balling and 
Sobieszczanski-Sobieski (1994). The CO strategy has been 
successfully applied to a number of different design problems 
including the design of a single-stage-to-orbit launch vehicle 
(Braun, et. al, 1996b). Tappeta and Renaud (1997) extended this 
approach and developed three different formulations to provide 
for multiobjective optimization of non-hicrarchic systems. 

Collaborative Optimization (CO) is a two level optimization 
method specifically created for large-scale distributed-analysis 
applications. The basic architecture of collaborative 
optimization (Braun, et. al, 1996a) is shown in Fig. 3. The CO 
framework facilitates concurrent design at the discipline design 
level. To achieve the concurrency in the subspace level, 
auxiliary design variables (x aux )jj arc introduced as additional 
design variables to replace the coupling variables so that the 
analyses in disciplines / and j can be executed concurrently. 
Compatibility constraints (or discrepancy functions) d are added 
to ensure consistency such that (x al(X ) i j=y ij . Compatibility 
constraints are usually in the form of equality constraints. 

The system level optimizer attempts to minimize a system 
level objective function F while satisfying all the compatibility 
constraints d. System level design variables x svs consist of not 
only the shared variables x sfl but also the auxiliary x aux 
variables. These variables arc specified by the system level 
optimizer and are sent down to subspaces as targets to be 
matched. Each subspace, as a local optimizer, operates on its 
own set of design variables x ssi with the goal of matching target 



Figure 3. Basic Collaborative Optimization 
Architecture (Braun et. al, 1996a) 


values posed by the system level as well as satisfying local 
constraints The matching can be attained by minimizing the 
discrepancy d l between some of the local design variables and/ 
or local states and their corresponding target values, in other 
words, the objective functions at subspace level are identical to 
the system level (compatibility) constraints. This formulation 
allows the use of post-optimal sensitivities calculated at the 
subspace optimum to be used as the gradients of the system 
level constraints. This important feature improves the overall 
efficiency of CO by eliminating the need to execute subspace 
analysis for the sole purpose of calculating system constraint 
gradients by finite differencing. 

DECISON-BASED COLLABORATIVE OPTIMIZATION 

Decision-Based CO (DBCO) Framework 

In the decision-based collaborative optimization (DBCO) 
framework developed in this research, the multidisciplinary 
enterprise model (Fig. 2) is decomposed along analysis 
boundaries into several subsystems. The method of 
collaborative optimization (CO) is used to determine the 
optimal design of this complicated model, where an optimizer is 
integrated within each subsystem or discipline. The resulting 
decision-based collaborative optimization (DBCO) framework 
rigorously simulates the existing relationship between business 
and engineering in multidisciplinary systems design as shown in 
Fig. 4. A brief discussion about this framework is given below. 

System Level Optimization 

In the decision-based collaborative optimization framework, 
the business decisions are made at the system level. As 
discussed before, the goal of decision-based design optimization 
is to maximize the expected utility of net revenue (E(u(NR))) of 
the engineering artifact being designed. The system level 
optimizer in the decision-based collaborative optimization 
framework attempts to increase expected utility of net revenue 
while satisfying compatibility constraints d. According to the 
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analyses in the business discipline and the subspace 
optimization results, the system level optimizer determines 
price P and establishes a set of performance targets for shared 
design variables and auxiliary design variables x mtXJ 
including demand q and total cost C -p Attributes a are treated as 
auxiliary design variables which influence demand. Note that 
demand q and total cost C T are also among the auxiliary 
variables. Also note that since the design variable price P is only 
needed in the business discipline, it is not a target for any 
subspace. 

The business discipline is operated on directly by the system 
level optimizer. It is not further decomposed into demand and 
utility of profit subspaces. The reasons why the business 
discipline is dealt with in this manner are as follows: First of all, 
there is no two-way coupling between the analysis in the 
demand subspace and the utility of profit subspace. The demand 
q is fed forward into the Profit subspace. Secondly, the analyses 
involved in this discipline are relatively simple and 
straightforward (Eqs. (1) and (2)). Consequently it is relatively 
easy to obtain sensitivity information with respect to the profit 
(or net revenue) or the utility of net revenue. In the cases where 
an analytical demand model and an analytical utility model arc 
supplied, the sensitivity information is readily acquired through 
the application of the chain rule. 

The system level optimization problem in its standard form is 
given in Eq. (3). 


Minimize: F = - E (u(NR)) 

w.r.t. x sys 

Subject to: d i = 0 / = 1,2, , « vv 

( X J < x < ( x ) 

K sys'min ~ svs ~ K sys'max 

o .00. 

x sys = Ka- X „ur) 

P > 0 (3) 

Note that maximizing the expected utility E(u(NR)) is 
equivalent to minimizing its negative value. Note also that net 
revenue (NR) is determined by the system level design variable 
x sys and price P. Since system design variables x svs are posed as 

targets to the subspaces, the term x° sys is used in Eq. (3) so that 

ihey can be distinguished from subspace design variables in the 

* 

subproblcm formulations. The term d i refers to the optimal 
value of the discrepancy function d i obtained by the 
subproblems. The formulation of d x is discussed in the next 
section. 

Subspacc Level Optimization 

The subspace optimizer seeks to satisfy the targets sent down 

by the system level optimizer and reports the discrepancy d { 

back to the system level. Meanwhile the subspace optimizers 
are subject to local design constraints gp In the field of 



Figure 4. Decision-Based Collaborative Optimization 
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engineering design, the design constraints normally guard 
against failure or restrict the search to preferred region of the 
design space. One example of failure-related constraints is to 
require that “the axial load in a beam not exceed its buckling 
load”. The statement that “the mass of the beam should be less 
than 7 kg” is an example of constraints based on preference. The 
use of constraints to restrict the search to preferred region of 
design space is not recommended in decision-based design 
approach. Firstly, to impose a preferred space, the engineer must 
decide and quantify what level of behavior is unacceptable or 
undesirable (i.e., not preferable). This is a matter of decision- 
making and by imposing constraints of preferences, the designer 
is removing some degrees of freedom in the design process. The 
resulting system optimization may fail to identify the design 
with the best trade-off, especially when this constraint is active 
or near active at the optimal solution. In the example of the 
beam, a beam of mass greater than 7 kg is said to be 
unacceptable. But it is possible that a beam of 7.1 kg can 
support a much higher load than a beam of 7kg. If the goal of 
the optimization is to find a light beam that can support a large 
load, the beam of 7.1kg might be a better design. Yet if a 
constraint is set to ensure the mass of the beam to be no greater 
than 7 kg, the optimizer will not locate the beam of 7.1kg, even 
though it may be more preferred. 

Upon closer examination, undesirable behaviors (i.e., non- 
preferred region) are often undesirable because such behaviors 
lead to a decrease in the demand of the product and/or an 
increase in the cost of the product, which in the end, results in a 
decrease of profit. In DBD the market place is used to determine 
preferred region of the design space through demand and cost 
models, and therefore constraints related to undesirability or 
preference are eliminated. Therefore the local constraints in the 
decision-based collaborative optimization (DBCO) framework 
tend to be those constraints that guard against system failure. 
Other traditional engineering constraints related to consumer 
preference are eliminated from the subspace optimization and 
are instead incorporated in the demand model and/or cost 
model. 

The subspace optimization problem for discipline l of Fig. 4 
in its standard form is given below: 

Minimize: 0 2 

wrt - x ss l d \ = (K/,)i~K/,)|) 


The subspace must satisfy local constraints g l while attempting 
to minimize discrepancies in system level targets. Note that 
attribute targets arc imposed in the second and the third terms of 
the discrepancy function. 

System Level Constraints Gradient 

The gradient of system level constraints plays an important 
role in forming search directions for the system level 
optimization. As mentioned earlier, one important feature of CO 
is that post-optimality sensitivity analysis from the converged 
subspace optimization problem can be used to provide system 
level derivatives for compatibility constraints (Kroo et. al, 

1 994). As a result, both computational expense and numerical 
error are reduced. This is possible because the system level 
design variables are treated as parameters (i.e., targets) in the 
subproblems. Note that for a certain discipline /, depending on 
rhe contributing analysis involved, not all the system level 
design variables x are necessarily posed as targets to be 

matched. It is possible that only a subset of x svs , referred to as 
=;* )., is sent down as subspace i targets. Generally all the 

subsets for n ss subspaces are not mutually exclusive, i.e., their 
intersections exist. The set of system level design variables x sys 
is the union of all the subsets (x svs ) . The gradient of system 

level constraint d t with respect to the subset (x svs ) of the 

system level design variables sent down as targets to discipline i 
is given below in Eq. (5). The gradient of system level constraint 

7 with respect to those system level design variables which arc 
not imposed as targets for discipline / is apparently zero. 

♦ 

= -2iu;,) f -<o,i (5 > 

The term (a' vj ), refers to the vector formed by the converged 
optimal values of local variables and states in discipline i at the 
end of subspace optimization. The elements in (jc vv ) j - arc the 
optimal counterparts of the system level targets (x n . v ). . 

Test Problem 


Subject to; 




+ X (>’u _ (- v nur)u) 


k 2 


0 

Al 


> 

VI 

1 — («v 1 ) 

' •'>» nun 

1 v f.w ' max 

- v ,vi = 

l) 


(4) 


A preliminary application of the decision-based collaborative 
optimization framework has been tested on an Aircraft Concept 
Sizing (ACS) problem. This problem was originally developed 
by the MDO research group at the University of Notre Dame 
■; Wujek, Renaud, Batill, et. al., 1996). It involves the preliminary 
sizing of a general aviation aircraft subject to certain 
performance constraints. The design variables in this problem 
ore comprised of variables relating to the geometry of the 
aircraft, propulsion and aerodynamic characteristics, and flight 
regime. Appropriate bounds are placed on all design variables. 
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The problem also includes a number of parameters which are 
fixed during the design process to represent constraints on 
mission requirements, available technologies, and aircraft class 
regulations. 

The original problem has ten design variables and five 
parameters. The design of the system is decomposed into six 
contributing analyses. This problem has been modified by 
Tappeta (1996) to fit the framework of multiobjective coupled 
M DO systems. It is further modified in this research to be 
suitable for the case of decision-based design (DBD). For 
comparison, the rest of this section gives a brief description of 
the modified ACS problem by Tappeta. It will be referred to as 
the ACS problem from hereon. The DBD version of the ACS 
problem will be discussed in the following sections. 

The Aircraft Concept Sizing (ACS) problem has three 
disciplines (see Fig. 5): aerodynamics (CA1), weight (CA2) and 
performance (CA3). It can be observed from the dependency 
diagram that the system has two feed-forwards and there are no 
feed-backs between disciplines. Table 1 lists the design 
variables and their bounds in the ACS problem. Table 2 lists the 
usage of design variables as inputs to each discipline. It can be 
seen that there are five shared design variables (x\ ~ and x 7 ). 
Table 3 lists the parameters and their values. Table 4 lists the 
state variables and their relations with each discipline. Clearly 
there are two coupled states (y 2 and j> 4 ). Table 5 contains all the 
relevant information for the ACS problem in the standard MDO 
standard notation introduced earlier. 


Table 1: List of Design Variables in ACS Problem 



Name (Unit) 

L 

U 

X\ 

aspect ratio of the wing 

5 

9 

x 2 

wing area (ft 2 ) 

100 

300 

*3 

fuselage length (ft) 

20 

30 

*4 

fuselage diameter (ft) 

4 

5 

*5 

density of air at cruise altitude (slugs/ft 3 ) 

.0017 

.002378 

x 6 

cruise speed (ft/sec) 

200 

300 

X 1 

fuel weight (lbs) 

100 

2000 


Table 2: Input Design Variables of Each Discipline in ACS 



CA1 (Aero.) 

CA2 (Weight) 

CA3 (Perf.) 

*l 

V 

V 


x 2 

V 

v 

V 

x 3 

V 

v 


x 4 

V 

v 


x 5 


V 


x 6 


V 


x 7 


V 

V 


^Shaded cells in the tabic indicate shared variables. 


Wetted Area 1 1 

1 J . 

*1~*4 o 

j Lift/Drag I 

Aero. J 


y 2 


I y* 

Total Weight t 1 

L j 

T y I 1 

Prp5 1 1 j Empty Weight I 

Weight J 


y 4 


y 3 


, ! «**»«<= i — y 5 

x 2 ,x 7 •:> i j 

*— — — 1 

^6~^8 J Stall Speed 1 — ► y 6 

Performance 1 — — — — ^ 


Figure 5. Aircraft Concept Sizing Problem 
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Table 3: List of Parameters in ACS Problem 


Table 5: Design Vectors and Functions for ACS Problem 



Name 

Description 

Value 

P\ 

Npass 

number of passengers 

2 

Pi 

N c „ 

number of engines 

1 

Pi 

w cn 

engine weight 

197 (lbs) 

Pa 

w 

pay 

payload weight 

398 (lbs) 

Ps 

N zult 

ultimate load factor 

5.7 

P6 

Eta 

propeller efficiency 

.85 

Pi 

c 

specific fuel consumption 

.4495 (Ibs/hr/hp) 

Ps 

c u 

maximum lift coeff. of the wing 

1.7 


Table 4: Lists of States in ACS Problem* 



Name (Unit) 

Output From 

Input To 

y\ 

total aircraft wetted area (ft 2 ) 

CA1 


yi 

max lift to drag ratio 

CA1 

CA3 

yi 

empty weight (lbs) 

CA2 


T4 

gross take-off weight (lbs) 

CA2 

CA3 

-^5 

aircraft range (miles) i 

CA3 


>'6 

stall speed (ft/sec) 

CA3 



*Shaded cells in the table indicate coupling states 


The objective in the ACS problem is to determine the least 
gross take-off weight within the bounded design space subject 
to two performance constraints. The first constraint is that the 
aircraft range must be no less than a prescribed requirement, 
and the second constraint is that the stall speed must be no 
greater than a specified maximum stall speed. In standard form, 
the optimization problem is given below: 


Minimize: 

F 

= Weight = y 4 

Subject to: 

S\ 

n 

i 

N 

O'" 

IV 

c 






= ] J a " se ^ >0 


y$ 


‘'stalls = 70 ft/scc 

Range rcii = 560 miles 



Vector or Function 

Variables or Content 


X 

[xj, x 2 , X 3 , x 4y v 5 ,x 6 , a* 7 ] 


x sh 

[x,,x 2> x 3 ,x 4 ,x 7 ] 


X au.x 

goals for y 2 , y 4 


0 

X syx 

f 0 O 0 0 0 o o. 

[A-,, x 2 ,x 3 ,x 4 , x-j y y 2 ,y 4 J 


F 

?=y 4 


system targets to 
be matched 

/ 0 x r o o o o o, 

( X sys ) 1 “ l * 1 * x 2' x y x 4> y 2 1 



empty set 


(*v/,)l 

[x,,x 2 ,x 3 ,x 4 ] 

CA 1 

( X aux ) l 

empty set 

x ss\ 

[x h X 2 ,Xj,X A ] 


S\ 

empty set 


analysis 

= CA l[x,,x 2 ,x 3 .x 4 ] 


optimal target 
values 

U„), = [Jr,,x 2 ,x 3> 4 .^] 


system targets to 
be matched 

y o . ,00000 0, 

( x sysh = U,.X 2 ,Xj.X 4 ,X 7 ,>' 4 ] 


x 2 

[ x 5> x (,} 


( x si,h 

fX|,X 2 ,X 3 ,X 4 ,X 7 ] 


< x aiah 

empty set 

CA2 

x ss2 

[X|,X 2 ,X 3 ,X 4 ,X 5 ,X 6 ,X 7 ] 


Si 

empty set 


analysis 

lv 3 .>' 4 ] = 

CA 2 l.V|, x 2 , x 3 . x 4 , x 5 , x 6 , ,r 7 | 


optimal target 
values 

( x v,h = Up v 2 ,x 3 ,x 4 .x 7 ,.v 4 | 


system targets to 
be matched 

“ 1*2’ x 7* y 2’ > 4 ] 


- v 3 

empty set 


( x shh 

U 2 . x j] 

CA3 

( X aux)l 

L)’2’ >’ 4 ] 

X sx 3 

[x 2 , x 7 , y 2 , >’4] 


Si 

[S\-Sl\ 


analysis 

I Fs' y*6 1 “ L /< 3 [ x 2 , * 7 , y 2 , > 4 1 


optimal target 
values 

(•4,) 3 = U'2’4, y\,y\\ 
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Demand Model and Cost Model 

The approach of decision-based design requires engineers to 
not only focus on the product performance but also life-cycle 
costs as well as demand and the profit obtained over the life- 
cycle of the product. Thus it is very important to construct a 
proper demand model and a proper cost model for the product. 
The authors are aware that the task to build such models is not 
an easy one, and engineers are generally not trained for this 
task. Since this research is concentrated on the optimization 
aspect of decision-based design, it is reasonable to assume that 
other discipline experts have developed such demand and cost 
models and made them available to the optimizer. 

In the case of the Aircraft Concept Sizing (ACS) problem, 
neither a demand model nor a cost model was available from the 
previous studies. In order to apply the decision-based 
collaborative optimization framework of Fig. 4, a demand 
model and a cost model have been developed. These models are 
built in a way such that they agree with industry trends for this 
specific class of aircraft and that they lead to reasonable 
optimization behavior. Although they are by no means 
complete, they serve fairly well as concept models for the 
application of DBCO at the current stage. Figure 6 illustrates 
the demand and cost models of the ACS problem. Only the 
annual demand, annual cost and annual profit are considered in 
the current research. 


Demand Model 

The first step in building the demand model is to identify the 
attributes that influence the demand q of this aircraft. The 
conventional optimization ACS problem (Eq. (6)) tries to 
minimize gross take-off weight (y 4 ) while satisfying two 
performance constraints, one on stall speed (v 6 ) and the other on 
aircraft range (y 5 ). Closer examination reveals that the objective 
function and two constraints, imposed in the original problem, 
are based on the designer’s estimate of customer preference for 
weight, stall speed and range. In the decision-based design 
approach, it is more appropriate to treat these quantities (take- 
off weighty, aircraft range y 5 , and stall speed y 6 ) as attributes of 
demand. Hence there are no performance constraints in the 
DBD version of the ACS problem, and the goal of the 
optimization is to maximize profit. 

It is also assumed that customers are interested in the cruise 
speed (.v 6 ) of the aircraft as well as how much room they would 
have on the airplane. A new state variable, fuselage volume (y 7 ), 
s introduced to reflect the concern for passenger room on the 
aircraft. In all, there are five attributes of demand in the DBD 
version of ACS problem: take-off weight y 4 , aircraft range y 5 , 
stall speed y 6 , fuselage volume y 7 , and cruise speed .v 6 . Note that 
demand is also influenced by price P. 

The demand model developed is a multiplicative model: 



11 



q = q(a , P) = q B 


FK 
U-i , 


' price 


where a = {y 4 , y 5 , y 6 , » x 6 } 

n a = 5 
q B = 1200 


(7) 


The term q B represents a baseline demand, which is set to 1200. 
The number of attributes is denoted by n a . The affect of each 
attribute on the final demand is reflected by the demand 
coefficient c r Similarly the term c pnce denotes the demand 
coefficient of price P. The final demand q is the product of all 
demand coefficients and the baseline demand q B . 

Demand coefficients for each attribute and price arc 
developed by financial analysts and marketing personnel within 
the business discipline and vary with time. For the purpose of 
this study, they are assumed fixed with respect to time and are 
given in Fig. 6. The curves in Fig. 6 plot the coefficient of 
demand on the ordinate and the corresponding attribute (or 
price) on the abscissa. 

1 . gross take-off weight (y 4 ) 

The lower the take-off weight, the higher the demand; but 
an aircraft with a very light weight is not desired. 

2. aircraft range (y 5 ) 

The longer the aircraft range, the higher the demand; but 
after the range reaches more than 600 miles, there is no sig- 
nificant increase in demand when range increases. 

This formulation is not unlike the original performance 
constraint g 2 (Eq. (6)), where the coefficient is set to 1 
when aircraft range equals to 560 miles. 

3. stall speed Cv6) 

The lower the stall speed, the higher the demand; but a 
near- zero stall speed is not necessary. 

This formulation is not unlike the original performance 
constraint gj (Eq. (6)), where the coefficient is set to 1 
when stall speed equals to 70 ft/sec. 

4. fuselage volume (y 7 ) 

The larger the fuselage volume, the higher the demand. 

5. cruise speed (x 6 ) 

The faster the cruise speed, the higher the demand. 

6. price ( P) 

The lower the price, the higher the demand; if the aircraft is 
sold for free (P= 0), the demand approaches infinity. 

Cost Model 

It is assumed that all costs (Cp) related to the production of 
the aircraft can be divided into two categories: fixed cost Cp and 
variable cost Cy (Krajewski and Ritzman, 1999). The fixed cost 


C F is the part of the total cost C T that remains consistent 
regardless of changes in the amount of product produced, for 
example, the annual cost of renting or buying equipment or 
facilities. The variable cost Cy is the portion of the total cost C T 
that varies directly with quantity of product produced, such as 
cost per unit for material and labor. If we assume the quantity of 
product produced and sold per year is equal to the demand q for 
the product per year, the total cost C^per year is: 

C r - C F + qC v (8) 

It is assumed that the variable cost Cy in the ACS problem is 
dependent on five of the seven design variables including wing 
area (x 2 ), fuselage length (x 3 ), fuselage diameter (.* 4 ), cruise 
speed (x 6 ) and fuel weight (x 7 ). A variable cost parameter (py) t 
is assigned to each cost-related variable to represent the portion 
of variable cost (per unit) associated with each variable. The 
total variable cost per unit is the sum of all variable cost 
parameters: 


n CU5f 

= I ( P y). (9) 

i = 1 

where n cos{ is the number of cost-related variables. The 
guideline for assigning variable cost parameters is: the larger the 
variable, the higher the cost. The curves in Fig. 6, associated 
with each cost-related variable, plot the variable cost parameter 
-unit: 10,000 dollars) on the ordinate and the corresponding 
cost-related variable on the abscissa. The step jumps in the 
curves represents the need to purchase (or rent) and/or install 
new equipment (or facilities) when the size of the aircraft 
exceeds existing production capabilities. 

Substituting Eq. (9) into (8), the model of the total cost Cp is: 


c r = Cp + q £'</>,.), U0) 

l - 1 

md the unit cost C u can be obtained by dividing both sides of 
Eq. (10) by demand q\ 


Q ^cosi 

C u = ~f + I (Pi’), (U) 

Note that the number of product produced (q) may have a 
discounting effect on the variable cost parameters py For 
instance, usually the cost per unit for material will decrease 
when the total amount of material bought increases. Thus a q- 
discounting option has been included in the determination of 
variable cost parameter in the cost model. 


12 



Note 

The demand and cost models developed in this paper are by 
no means complete. They are conceptual and rather simplistic. 
Future work on the modification of these models will include 
(and not limit to) the following issues: 

1 . gross take-off weight (y 4 ) 

It has been pointed out to the authors that to a customer 
higher gross weight is actually desirable because it leads to 
longer aircraft range. Meanwhile higher gross weight leads 
to higher manufacture cost. Therefore a modified demand 
model would include the gross take-off weight as a slightly 
favorable feature (i.e., the higher the take-off weight, the 
higher the demand). On the other hand a modified cost 
model may include the gross take-off weight as a strongly 
negative factor (i.e., the higher the take-off weight, the 
higher the cost). 

2. aspect ratio of the wing (X[) 

Increasing aspect ratio will increase the wing structural 
weight, which will in turn lead to an increase in the aircraft 
gross take-off weight, thus add to the total cost. A modified 
cost model would include aspect ratio as another negative 
factor. 

3. price (P) 

It has been brought to the authors’ attention that in the real 
world, due to the maintenance requirements such as insur- 
ance and hanger, the demand will not approach infinity if 
the aircraft is given out for free. A modified demand model 
will address this issue by assigning a definite number to the 
demand coefficient for price when price is set to zero. This 
definite number will be associated with the maximum 
demand possible for this aircraft. 

DBCO Formulation 

The decision-based collaborative optimization framework has 
been applied to the DBD version of the Aircraft Concept Sizing 
(ACS) problem. This application is a preliminary study, and 
focuses on the collaborative optimization feature of the DBCO 
framework. The issues of propagated uncertainty arc neglected 
in this study. The utility of profit is assumed to be the profit 
itself. Hence the objective of the resulting deterministic 
optimization is to maximize profit (or net revenue). During the 
optimization, the demand q is treated as a continuous variable, 
rather than an integer. At the end of the system optimization, q 
is rounded to the nearest integer. 

Two additional disciplines are added in the DBD version of 
the ACS problem: cost (CAc) and business (CAb). Price P is a 
new design variable and a new state variable (fuselage volume 
y 7 ) is introduced. Table 6 provides the list of input design 
variables to each discipline in the DBD version of ACS 
problem. Clearly design variable x 6 (cruise speed) enters the set 
of shared variables. Table 7 lists the states y , demand q , total 


cost C T , and net revenue NR. It also depicts how they are related 
:o each discipline. The set of coupling variables expands to 
include five additional members: V5 (aircraft range), v 6 (stall 
>peed), >’7 (fuselage volume), q (demand) and Cj (total cost), 
fable 8 contains all the design vector information for the DBD 
version of subspaces l (aerodynamics), 2 (weight), 3 
performance) and c (cost) of the ACS problem in MDO 
standard notation. Table 9 contains the design vector 
information of subspace b (business) since it is operated directly 
by the system level optimizer. The difference between the DBD 
version of ACS problem and the modified ACS problem by 
Tappeta (1996) can be clearly observed by comparing Tables 6, 
7, 8 and 9 with Tables 2, 4 and 5, respectively. 


Table 6: Input Design Variables to Each Discipline in ACS 
(DBD version)* 



CA1 

(Aero.) 

CA2 

(Weight) 

CA3 

(Perf.) 

CAc 

(Cost) 

CAb 

(Busin.) 

x \ 

4 

4 




\ x 2 

4 

4 

4 

4 


*3 

4 

4 


4 


X A 

4 

4 


4 


x 5 


4 




x 6 


4 


4 

4 

X 1 


4 

4 

4 


P 





4 


'Shaded cells indicate shared variables. 


Table 7: Lists of States in ACS Problem (DBD version)* 



Output From 

Input To 

y 1 

CA1 


T2 

CA1 

CA3 

V 3 

CA2 


>’4 

CA2 

CA3, CAb 

>’5 

CA3 

CAb 

>'6 

CA3 

CAb 

>’7 

CAI 

CAb 

<7 

CAb 

CAc 

C r 

CAc 

CAb 

NR 

CAb 



'Shaded calls indicate coupling states. 
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Table 8: Design Vectors for ACS Problem 
(DBD version) 



Vector or Function 

Variables or Content 


.V 

[.t], X 2 , X 3 , X 4> X5, X5, X-j, P] 


*sh 

[x,,.v 2 ,x 3 ,x 4> x 6 ,x 7 ] 


A aux 

goals for >>2, y 4 , y 5 ,y 6 , y 2 , C r , q 


o 

*,y. 

, O O 0 0 0 0 

[x t ,x 2 ,x } , x 4 ,x 6 ,x 7 , 

0 O 0 O O ~,0 0 , 

y^y^y^y^ y^c T ,q ] 


F 

F=-NR 

CA 1 

system targets to 
be matched 

.O. -OOOOOO, 

= 1*1 , X 2 , Xy x 4 , y 2 , y 7 ] 


empty set 

( x sh ) 1 

[x,,x 2 ,x 3 ,x 4 ] 

(*auxh 

empty set 

x ss 1 

[A|,X 2) X 3 ,X 4 ] 

8 \ 

empty set 

analysis 

[y v y 2 >yi\ = oiiix,,x 2 ,x 3 ,x 4 ] 

optimal target 
values 

= [x,,x 2 , x 3 , x 4 ,y 2 ,y 7 ) 

CA 2 

system targets to 
be matched 

.0. , 0 0 0 0 0 0 o, 

(X,y ,) 2 l x v X 2' x y x 4 , x b , x 7 , v - 4 1 

*2 

[x 5 ) 

(*v/i)2 

[x h x 2 ,xyx 4 ,x b ,x 7 ] 

(• X anxh 

empty set 

x s.s 2 

[*| , x 2 , x 3> x 4 , x 5 , x 6 , v 7 ] 

82 

empty set 

analysis 

[>' 3 . V 4 I = 

CA 2 [x,, x 2 , x y .v 4 , v 5 , v 6 . a 7 ] 

optimal target 
values 

( x ssh = x 3 , x 4 , x 4 , .r 7 0' 4 ] 

CA 3 

system targets to 
be matched 

. 0. ,00000 0, 
(x vvv ) 3 = [x l ,x v y 2 ,y 4 ,y s ,y b ] 

a '3 

empty set 

3 

[x 2 , x 7 ] 


L>’2. V4] 

A v.v 3 

{x 2 ,x 7 ,y 2 ,y 4 \ 

8 l 

empty set 

analysis 

l^ 5 . ,V 6 1 = CA 3 [x 2 , v 7 , v 2 , ,v 4 ] 

optimal target 
values 

(Xssh = [x 2 , x 7 . y 2 . y 4 , y 5 , y 6 1 


Table 8: Design Vectors for ACS Problem 
(DBD version) (Continued) 


system targets to 
be matched 

,o x ,0 0 0 0 0 0^,0, 

= [x 2 . x 3 . x 4 . x 6 . x 7 , q ,C T ] 

v c 

empty set 

(" V aV;)c 

[x 2 ,x 3 ,.v 4 ,x 6 ,x 7 ] 

( A auOc 

w 

A vvc 

[x 2 , x 3 , x 4 , x 6 , x 7 . q] 

8c 

empty set 

analysis 

C T = CAc[x v x y x 4 , x 6 , .v 7 , q) 

optimal target 
values 

(Xssh = [XyXyX 4 ,X b .X 7 ,q ,C r I 


Table 9: Design Vectors in the subspace b (business) 
for ACS Problem (DBD version) 


targets sent down 
to other subspaces 

, O x r 0 0 0 O 0 0 ^ 0 , 

(Xsysh = l r 6- - v 4’ y 5 ’ - V 6- >'?- 9 • C T 1 

A b 

[P] 

( A j/i )b 

[x 6 ] 

( A auv)b 

L>' 4 ’ >' 5 - >’6> >'7» C T\ 

'Vv.vb 

r O 0 0 0 0 0 r, , 

[-v 6 . y 4 , v s , y 6 , y 7 , C T ,P] 

8b 

empty set 

analysis 

lq°,NR) = 

^ a i r o 0 0 0 O /-i 0 jj i 

CAb[x b , y 4 , v 5 , y 6 , y v C T , P\ 


The system level optimization problem, for this application, 
in its standard form is detailed in Eq. (12). 


System Level Optimization 


Minimize: 

, o 

w.r.t. jr 


F - -NR 


Subject to: 


where 


d\ = 0 
d\ = 0 
d\ = 0 
/. = 0 

{ V° ) < JC° < ( v° ) 

v xy. r J nun ~ sys — v svs 'max 

o r o o o o o o 

x sys = [ Xi'X 2 ,XyX A ,X 6 ,X 7 . 

O O O O 0 ^ O , 

y f 2 ’ >4’ r 5 - r 6 > » c r ^ 

<7 = 40 * 4 ’ >' 5 ’ > ! 6- >* 7 > A 6’ 

NR = P q°-Cj 


( 12 ) 
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Note that the system level optimizer calls the business discipline 

directly to obtain demand q° and the system level objective NR. 
There are thirteen system level design variables and four 
compatibility constraints that arc evaluated in subspace 1, 2, 3 
and c. 

The subspace optimization problems for each discipline in 
their standard forms are given by Eqs. (13) ~ (16). 

Subs pace 1 (Aerodynamics) Optimization 

Minimize: 0 2 o 2 

w.r.t. x sxX </, — ( -V , — .Y j ) +(x 2 -x 2 ) 


o 2 o 2 

+ (* 3 -* 3 ) +(x 4 -jc 4 ) 


+ (y 2 - v 2 ) +(y 7 -v 7 ) 

Subject to: (x ,) <x . < (x ,) 

j v A.n 'nun •' s i ' aa i 'max 

where x xxl = [x x ,x 2l x y x A ] 

[ y 1 1 v 2 > y 7 ) = C/t l [,V|, a'2» ^ 3 , a' 4 ] 

Subspace 2 (Weight) Optimization 


Minimize: 
w.r.t. x xx2 


d 2 = (.Vj -.v°) + (* 2 -jr 2 ) 


(13) 


Subject to: 
where 


+ (x 3 -x 2 ) + (x 4 -x 4 ) 

o 2 o 2 

+ (>4-r4) 2 

x .is2 = [*,, * 2 . * 3 . * 4 . * 5 ' * 6 - * 7 ! 

[y 3 , y 4 J = C/t2[*,, * 2 , * 3 , * 4 , x 5 , x b , x 7 \ 


Subspace 3 (Performance) Optimization 


Minimize: 
w.r.t *„ 3 


dy (**2 — *2 ) + (* 7 — * 7 ) 


(14) 


2 2 
+ (W 2 -JS) + (> 4 " >4) 

^(>5-V5r+0* 6 -V6) 2 
Subject to: U„ 3 < -v, v3 < (-v v ,. 3 ) m(it 

where x ss3 = [x 2 , .v 7 , y 2 , v 4 ] 

1 J’5> J7> 1 = C-4 3 l v 2’ a 7’ v 2’ ^4 1 


(15) 


Subspace c (Cost) Optimization 

Minimize: 0 2 0 2 

wrtx sxc d c = (x 2 -x 2 ) +(x 3 -x 3 ) 

O 2 0 2 

+ (*4 -*4) + (* 6 -* 6 ) 
o 2 

+ (*7 -*7> 

+ (<?-<7°) +(C r -C r °) 
Subject to: S *„ c < 

where x xxc = [x 2 , * 3 , * 4 , x 6 , * 7 , q } 

C T = C4c[* 2 , .v 3 , * 4 , .r 6 , x 7 , <7] 


Mote that other than variable bounds, there are no local 
constraints for the subspace optimization problems. 


Optimization Results & Discussion 

A Sequential Quadratic Programming (SQP) method was 
ised for optimization in both the system level and the subspace 
optimization. The SQP solver, fmincon, was obtained from the 
Vf atlab Optimization Toolbox. The program converged to an 
optimum after thirty-seven system level iterations. The optimal 
solution is listed in Table 10. Figure 7 shows the system level 
optimization history of convergence of the system level 
objective function (negative of profit, in subplot 8), the 
convergence history of the four compatibility (discrepancy) 
♦ • * * 

constraints (d |# d 2 , d 3 , and d c , in subplots 9~12), and the 
convergence history of the seven system level design variables 
t cruise speed x° 6 , aircraft range y° , stall speed y° e , fuselage 

volume y 7 , price P, demand q° and total cost C^, in subplots 

1 ~7). The abscissa of each subplot is the number of system level 
iterations. Note that the value of profit (not the negative of 
profit) was plotted in subplot 8 for easy reading. For the same 

reason unit cost C° was plotted instead of total cost C 7 ° in 
subplot 7. 

As can be seen from the convergence plots, the system level 
-optimizer tries to minimize both the negative of profit and the 
constraint violations simultaneously. At the beginning of the 
optimization, the system level optimizer sets targets high for 
price, high for the levels of performance (to ensure high 
demand) and low for cost based on the results of the business 
malyses. However these targets conllict with one another and 
lead to a large discrepancy at the subspace level. Thus the 
system level optimizer, while trying to keep profit as high as 
possible, was forced to lower price, downgrade performance 
and tolerate higher cost so that the subspace discrepancy could 
be reduced. Gradually the system level optimizer found the best 
;rade-off among the targets and reached a consistent optimal 
design. The optimization history observed in the ACS problem 



resembles the existing relationship between business and 
engineering in multidisciplinary systems design. 

The demand model and the cost model play an important role 
in the decision-based design approach. In order to illustrate the 
influence of the demand and cost models, a conventional all-at- 
once optimization was performed according to the problem 
formulation in Eq. ( 6 ). The conventional optimum obtained is 
also listed in Table 10. Note that fuselage volume (j ; 7 ) at the 
conventional optimum is determined by the optimal fuselage 
length (JC 3 ) and optimal fuselage diameter (jc 4 ). It can be 
observed that the conventional optimum outperforms the DBCO 
optimal design on lower weight (v 3 , y 4 ). However it possess 
poor characteristics in many aspects such as smaller aircraft 
range ( 75 ), higher stall speed (y 6 ) and smaller fuselage volume 
(y 7 ). Such an outcome is no surprise since the main concern of 
the conventional ACS problem is to minimize take-ofT weight, 
while the DBD approach takes into account other performance 
attributes, because of the DBD objective of maximizing profit. 

If we assume that the aircraft configuration at the 
conventional optimum design will be sold at the same price as 
the DBD optimum design, the demand, cost and profit of the 
conventional product can be obtained according to the demand 
model (Eq. (7)) and cost model (Eqs. (10) and (11)) developed 
earlier. These values are listed in Table 10 in parentheses 
because of the assumption. Notice that the unit cost of 
conventional optimal design is lower than the unit cost of the 
DBD optimal design. However the poor performance attributes 
cause the demand for conventional optimal design to be much 
lower than the DBD optimal design. Hence the DBD optimal 
design leads to higher profit. 


Table 10: Optimal Solutions for ACS Problem 



Name (Unit) 

DV 

Bounds 

DBCO 

Optimum 

Convcn. 

Optimum 

*1 

aspect ratio of the wing 

5-9 

7.968 

5 

*2 

wing area (ft 2 ) 

100-300 

230.3 

176.53 

- v 3 

fuselage length (ft) 

20-30 

21.927 

20 

*4 

fuselage diameter (ft) 

4-5 

4.1871 

4 

X S 

density of air at cruise 
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Figure 7. System Level Convergence Plots for ACS Problem (DBD version) 
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CONCLUSIONS 

In this research a Decision-Based Collaborative Optimization 
(DBCO) framework which incorporates the concepts of 
normative decision-based design (DBD) and the strategies of 
Collaborative Optimization (CO) has been developed. This bi- 
level non-deterministic optimization framework more 
accurately captures the existing relationship between business 
and engineering in multidisciplinary systems design. The 
business decisions are made at the system level, which result in 
a set of engineering performance targets that disciplinary 
engineering design teams seek to satisfy as part of subspace 
optimizations. The objective of the Decision-Based 
Collaborative Optimization (DBCO) is to maximize the 
expected von Neuman-Morgenstcm (vN-M) utility of the profit 
or net revenue (NR) of a product. 

A preliminary application of this approach (deterministic 
case) has been conducted on a multidisciplinary test problem 
named the Aircraft Concept Sizing (ACS) test problem. 
Conceptual demand and cost models have been developed. The 
corresponding optimization results are discussed and compared 
with the conventional optimization solutions. 

Future work is being targeted towards a non-deterministic 
implementation of the DBCO framework in which the issues of 
propagated uncertainty in such a bi-levcl optimization 
framework will be addressed (Gu and Renaud, 2001). A 
conceptual utility model for the net revenue (NR) will be 
adapted from the literature in the field of decision-based design 
(DBD). The conceptual demand and cost models developed in 
this paper will be modified to better reflect the real world. 
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